TITLE OF THE INVENTION 
DATA CONVERSION APPARATUS AND METHOD 



BACKGROUND OF THE INVENTION 

FIELD OF THE INVENTION 

The present invention relates to a data 
conversion apparatus and its method and, more 
particularly, ' to a data conversion method suitable for 
digital arithmetic operations, and a data conversion 
apparatus using the method. 

DESCRIPTION OF THE RELATED ART 

The volume of image data digitally processed by 
computers and the like is increasing each year. Also, 
higher quality upon forming a color image is 
increasingly required each year. In order to form a 
high-quality color image, color conversion is 
indispensable, and requires high-speed processing in 
addition to high precision and implementation of 
flexible conversion characteristics. 

Since color information generally forms a 
three-dimensional space, color conversion determines 
the correspondence of color information from a given 
color space to another color space. Many schemes for 
attaining such conversion are available. Among these 
schemes, color conversion that combines a look-up table 
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(LUT) and interpolation (Japanese Patent Laid-Open 
Nos. 53-123201 and 8-237497) is prevalently used. Also, 
color conversion that combines an LUT and interpolation 
includes various schemes. In consideration of the 
required data size, computation volume, continuity of 
outputs among unit rectangular hexahedra, gray line 
interpolation characteristics, and the like, 
tetrahedral interpolation disclosed in Japanese Patent 
Laid-open No. 53-123201 is most suitable. 

However, the tetrahedral interpolation disclosed 
in Japanese Patent Laid-Open No. 53-123201 can be 
directly applied only when the unit rectangular 
hexahedron is a regular hexahedron. Upon implementing 
color conversion by a digital computation processing 
apparatus, the limitation of a unit rectangular 
hexahedron to a regular hexahedron cannot be ignored in 
implementation of a conversion processing ap^^aratus. 

This limitation will be described in detail below, 
in a conversion processing apparatus that performs 
digital processing, in order to convert all unit 
rectangular hexahedra into regular hexahedra, the grid 
spacing is limited, and the conversion precision and 
data size (the number of grid points) cannot be 
optimized. For example, if 'input data is 8-bit data (0 
, to 255),' the grid spacing must be set at one of 85 (the 
number of grid points = 4), 51 (6),- 17 (16), 15 (18), 5 
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(32), and 3 (86) to convert all unit rectangular 
hexahedra into regular hexahedra. 

On the other hand, since the color space is not 
uniform in the entire area, it is effective for 
improving the conversion precision of a specific area 
to intentionally set a small grid spacing of that area. 
For example, upon conversion from RGB to CMYK, a gray 
level drop between grids arising from an undercolor 
removal (UCR) process readily occurs in a dark area, 
i.e., (R, G, B) - (0, 0, 0). However, this problem can 
be effectively solved by setting a small grid spacing 
of that area. Japanese Patent Laid-Open Nos. 7-131668, 
10-70669, and the like disclose schemes for improving 
the conversion precision of a specific color area by 
setting a small grid spacing. However, the schemes 
described in these references require an area 
discrimination process for changing the processes 
depending on the areas of input color information, and 
are inferior to the technique disclosed in Japanese 
Patent Laid-Open No. 53-123201 above in terms of the 
processing speed and circuit scale. 

SUMMARY OF THE INVENTION 
It is an object of the present invention to 
provide a data conversion method suitable for digital 
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computation processing, and a data conversion apparatus 

using the method. 

In order to achieve the above object, the 
preferred embodiment of the present invention discloses 
a data conversion method comprising the steps of: 
outputting a value which represents distance from an 
input value to a grid point of a look-up table, and is 
normalized by a sufficiently large value, using the 
look-up table; and executing data conversion of the 
input value by interpolating the value obtained by the 

look-up table. 

It is another object of the present invention to 
provide a recording medium used in the data; conversion. 

In order to achieve the above object, the 
preferred embodiment of the present invention discloses 
a computer readable medium recorded data which is used 
in a data conversion process, the data comprising: 
table data for outputting a value which represents 
distance from a grid point of a look-up table to an 
input value, and is normalized by a sufficiently large 
value, with respect to the input value; and data 
representing a computation for executing data 
conversion of the input value by interpolating the 
value obtained by the look-up table, using the value 
25 obtained by the table data. 
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It is still another object of the present 
invention to provide an image processing apparatus and 
method that use the data conversion. 

In order to achieve the above object, the 
5 preferred embodiment of the present invention discloses 
an image processing method comprising the steps of: 
selecting a plurality of grid points on the. basis of 
input data; obtaining values, which represent distances 
between the selected grid points and the input data, 
10 and are normalized by a predetermined value; and 

executing interpolation on the basis of the obtained 
values and data of the plurality of grid points. 

Other features and advantages of the present 
invention will be apparent from the following 
15 description taken in conjunction with the accompanying 
drawings, in which like reference characters designate 
the same or similar parts throughout the figures thereof. 

BRIEF DESCRIPTION OF THE DRAWINGS 
20 Fig. 1 shows the entire two-dimensional LUT; 

Fig. 2 shows a given unit rectangle extracted 
from the two-dimensional LUT shown in Fig. 1; 

Fig. 3 is a view showing the principle of linear 
interpolation in case of one dimension; 
25 Fig. 4 is a view for explaining the contents of 

an X-u' table; 
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Fig. 5 is a view for explaining the relationship 
between x and distance to a grid; 

Fig. 6 shows the overall image of a 
three-dimensional LUT; 
5 Fig. 7 shows a unit rectangular hexahedron as a 

part of the three-dimensional LUT shown in Fig. 6; 

Fig. 8 is a block diagram showing the arrangement 
of a conventional apparatus; and 

Fig. 9 is a flow chart showing a conversion 
10 process. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 
The preferred embodiment of a data conversion 
apparatus and method according to the present invention 

15 will be described in detail hereinafter with reference 
to the accompanying drawings. 
[Problems in Tetrahedral Interpolation] 

In order to clarify problems posed when the 
interpolation (tetrahedral interpolation) disclosed in 

20 Japanese Patent Laid-Open No. 53-123201 is implemented 
by a digital computation processing apparatus and is 
applied to an LUT in which the grid spacings are not 
equal in all dimensions, the following explanation will 
be given taking two-dimensional interpolation as an 

25 example. As a scheme used in the interpolation to be 
described below, tetrahedral interpolation in a 
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three-dimensional space disclosed in Japanese Patent 
Laid-open No, 53-123201 is generalized and adapted to a 
two-dimensional space, and problems in the following 
example become more conspicuous in three-dimensional 
5 tetrahedral interpolation- Note that the description 

will be made in reference to a two-dimensional space to 
simplify the drawings, and to make it easy to 
understand the drawings. 

Fig. 1 shows the entire two-dimensional LUT, and 

10 Fig. 2 shows a given unit rectangle extracted from the 
two-dimensional LUT. Note that a unit rectangular 
hexahedron in three-dimensional tetrahedral 
interpolation becomes, a unit rectangle in the 
two-dimensional space . 

15 Interpolation in the two-dimensional space is 

done by checking which of triangles POO-PlO-Pll and 
POO-POl-Pll includes an input point, and . multiplying a 
value (grid value) at the vertex of the triangle that 
includes the input point by a weight. Assuming that 

20 POO - P(XO, YO), PIO = P{X1, YO), POl = P(XO, Yl) , and 
Pll = P(X1, Yl) (XI - XO Yl - YO), a process for 
calculating an output P for an input point (X, Y) that 
satisfies XO < X < XI and YO < Y < Yl will be explained 
below, 

25 If u and v are defined by: 




...da) 

(XI -xo) 

v = -fczYO) ...(lb) 
(Yl-YO) 

which of the triangles includes the input point (X, Y) 
is determined by comparing u and v, and a grid value 
used in interpolation differs. 

Case 1: if u > v, the input point (X, Y) is 
included in triangle POO-PlO-Pll. Hence, the output P 
is given by: 



(XI - x){Yi - YO)poo + {{X - xo)(n - YO) - (XI - xo)(y - yo)}Pio + (xi - xo){y - yo)pi i 
^ffff (xi-xo)(n-KO) 

10 ...(2a) 

Case 2: if V > u, the input point (X, Y) is 
included in triangle POO-POl-Pll- Hence, the output P 
is given by: 

~xi - XQ-jn -r)m+ {{xi -xo)l.r-YO)-ix- xojn - ro)]m t (x - xoxn - ro)f 1 1 



^~ (xi-xo)(n-yo) 




in line segment POO-Pll. Hence, the output P is given 
by : 

(XI - X)POO + (X - X0)P1 1 ...(2c) 



^ XI- XO 



20 Since equation (2c) is equivalent to equation (2a) or 
(2b) when u = v, case 3 can be included in case 1 or 2 
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If (XI - XO) = (Yl - YO), i.e., if the unit rectangle 

is a square, equations {2a) and (2b) are greatly 

simplified as: 

(XI - X)POO + {(X - XO)- (Y - Y0)>>1 0 + (Y - Y0)P1 1 

Xl-XO ... (2a') 



5 p _ (Yl - Y)POO + {(Y - YO)- (X - X0);^1 0 + ((X - X0)P1 1 

Xl-XO ...(2b') 

In this manner, when the unit rectangle is a 
square in tetrahedral interpolation in a 
two-dimensional space, calculations of the output P do 
not require any multiplication like (XI - XO) (Yl - Y) , 

10 the computation volume can be greatly reduced. 

Likewise, even in tetrahedral interpolation in a 
three-dimensional space, if the unit rectangular 
hexahedron is a regular hexahedron, the computation 
volume upon calculating the output P can be reduced. 

15 However, if the unit rectangular hexahedron , is not a 
regular hexahedron, calculations of one output P 
require multiplications like (XI - XO) (Yl - YO) (Zl 
- ZO) . Hence, the computation volume increases 
considerably, and an apparatus that can attain 

20 high-speed processing and has a large circuit scale is 
required. 

Of course, if floating point computations are 
made in digital processing, equation (2a) is rewritten 
as : 
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P = (l-XX)P00-h(XX-YY>10 + YYxPll . . . (2a") 

for 

XX is a real number given by XX = (X - X0)/(X1 

- XO) (0 < XX < 1) ... (3a) 
5 YY is a real number given by YY = (Y - Y0)/(Y1 

- YO) (0 < YY < 1) ... (3b) 
Hence, even when (XI - XO) ^ (Yl - YO) , equation (2a) 
can be simplified as equation (2a ') . 

However, floating point computations require a 

10 considerably longer processing time than integer 

computations, and also a huge circuit scale. For this 
reason, floating point computations of equation (2a") 
often increase the processing time and circuit scale 
compared to integer computations of equation (2a) , and 

15 high-speed interpolation cannot be implemented by a 
simple circuit. 
[Linear Interpolation] 

The basic principle of this embodiment will be 
explained below taking as an example a case wherein 

20 this embodiment is applied to simplest, linear 
interpolation . 

Given PO = P(XO) and PI = P(X1), as shown in ■ 
Fig. 3, if linear interpolation is done, an output P 
for an input point X that satisfies XO < X < XI is 

25 given by: 
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(X1-X)P0 + (X-X0)P1 

fUl ' XI -xo 



However, in equation (4), individual values (XI 
- XO), (Xl - X), and (X - XO) have no bearing, but the 
ratios of two values, i.e., (XI - X) / (XI - XO) and (X 
5 - X0)/(X1 - XO) have bearing as weights. Hence, if u 
is calculated using an arbitrary value L by: 

^ ^ x-xo , 



u = 



(5) 



XI -XO 

— ■ equation (4) is rewritten using u as: 
^^ (L-u)PO-mxPl 

10 When equations (5) and (6) are calculated by 

floating point computations, u becomes a real number, 
and equations (6) and (4) are equivalent to each other. 
However, it is inappropriate to make floating point 
computations,' as described above. For this reason, 
15 when equivalent processes are done by integer 
computations, u obtained by equation (5) is an 
approximate value. Hence, computations done by an 
actual apparatus are accurately expressed by: 

a = ROUND( ^'^^^ l] ... (5') 

Ul-XO J 
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where ROUND (x) is a function of rounding off after the 
decimal point of real number x. 

Hence, the output P obtained by equation (6') 
does not perfectly match that obtained by equation (4). 
5 However, if a value sufficiently larger than {XI - XO) 
is set as L, that error is negligible in practice. In 
addition, since L can be independently set except for 
(XI - XO) and its value, if a power of 2 is set as L, a 
division in equation (6') can be implemented by a 
10 right-shift computation. In general, the shift 

computation can be processed by a simpler circuit and 
at higher speed than the division. Therefore, by 
setting a power of 2 as L, reductions of the processing 
time and circuit scale can be realized. 
15 Since u' depends only on X, a table that outputs 

u' using X as an address, i.e., an X-u* table, is 
created prior to interpolation in actual processing. 
Fig. 4. shows the contents of the X-u' table as a graph. 
For the sake of comparison. Fig. 5 shows as a graph the 
20 relationship between X and distance to a grid, when L 
is not used. This X-u' table may be created 
immediately before execution of interpolation. If grid 
positions are determined in advance, the X-u' table may 
be stored in the apparatus in advance in place of the 
25 grid positions. 
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In this manner, if the X-u' table is prepared, 
the output P for the input X can be calculated by very 
simple computations using equation (6'). 
[Two-dimensional Interpolation] 

In linear interpolation, a division can be 
replaced by a shift computation, and larger effects can 
be expected when tetrahedral interpolation is used in a 
space not less than two dimensions . To demonstrate 
such effects, the following explanation will be given 
taking two-dimensional tetrahedral interpolation shown 
in Fig. 1 as an example. As a scheme used in the 
following description of the interpolation, tetrahedral 
interpolation in a three-dimensional space is applied 
to a two-dimensional space, and the effects in the 
following example are more conspicuous in the 
three-dimensional tetrahedral interpolation . 

As in u' given by equation (5) in case of one 
dimension, u' and v' in a two-dimensional space are 
given by: 



The relationships among u and v, and u' and v' are as 
follows : 



u'= INT 



X- XO ^ 

XI - XO 



(7a) 




(7b) 



r 

m 



if u > V, since uL > vL, INT(uL) > INT(vL), i.e., 
u ' > V ' ; 

if V > u, since vL > u-L, INT(vL) > INT(uL), i.e., 
V ' > u * ; and 

5 if u = V, since uL = vL, INT(uL) = INT(vL), i.e., 

u' = v' . 

As can be seen from the above description, by 
converting u and v into u' and v', the same 
relationship between u' and v' as that between u and v 

10 can be maintained although inequality sign may change 
to equality sign. Hence, since the triangle that 
includes the input point remains the same, it is 
appropriate to calculate the output P based on the 
relationship between u' and v' . Applying equations 

15 (7a) and (7b) to equations (2a) and (2b) yields: 
When u' > V, p, (L-u'>-00^(u^v»OH-VPn 

lOi^ When V > u.. p, (L-v»00^(V-u'>'01^u'Pll 



As shown in equations (8a) and (8b), 
interpolation for arbitrary grid spacings can be 
20 implemented by a computation volume equivalent to that, 
required when the grid spacings are equal to each other 
in equations (2a') and (2b'), i.e., when the unit 
rectangle is a square. Of course, when L is a power of 
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2, divisions by L in equations (8a) and (8b) can be 
implemented by right-shift computations. 
[Three-dimensional Interpolation] 

The processing sequence upon application of this 
embodiment to tetrahedral interpolation using a 
three-dimensional LUT will be explained below. Fig. 6 
shows the entire image of a three-dimensional LUT, and 
Fig. 7 shows a unit rectangular hexahedron as a part of 
the three-dimensional LUT shown in Fig. 6. 

Assuming that POOO = P(XO, YO, ZO), POOl = P(XO, 
YO, ZD, POlO = P(XO, Yl, ZO), POll = P(XO, Yl, Zl), 
PlOO = P(X1, YO, ZO), PlOl = P(X1, YO, Zl), PllO - P(X1, 
Yl, ZO), and Pill = P(X1, Yl, Zl), the effects obtained 
upon applying this embodiment to three-dimensional 
tetrahedral interpolation will be explained taking as 
an example the process for calculating an output P for 
an input (X, Y, Z) that satisfies XO < X < XI, YO < Y < 
Yl, and ZO < Z < Zl. 



u'= INt( ii_ 



xo ^ 

XO ] 



lip 



v'=INT 



Y-YO . 



Yl-YO 



w'= INT 



Z- ZO ^ 
Zl - ZD 



(9a) 



. (9b) 



(9c) 



As in the case of the two-dimensional 
interpolation mentioned above, since u' depends only on 
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X, an X-u' table as a table for obtaining u' using X as 
an address is created prior to interpolation. Likewise, 
a Y-v' table that obtains using Y as an address, and 
a Z-w* table that obtains w' using Z as an address are 
5 created prior to interpolation. These tables may be 
created immediately before execution of interpolation. 
When X, Y, and Z grid positions are determined in 
advance, the X-u ' , Y-v', and Z-w* tables may be stored 
in the apparatus in place of the grid positions. If 
10 these tables are created, the output P for the input (X, 
Y, Z) can be calculated very simply by: 

Whenii>v'>w' ,P = ^ '- ^ '- '- . . . (10a) 



4 



Whenu'>w'>v',P = ^ '- =^ ^- . . (10b) 

Li 

1,71. . . D - W )P000 + (w- 001 + ju' -^/)Pm + V' PI 1 1 ^ 
W/ienH'>M>v, P = -^^ ■ — . . . (10c) 



15 Wto»/>..^.f.''-''"^"'"'-''"^'^'-'-'^'""'^'^''"' ... (lOd) 



L 



. , n (L - V' )P000 + (v' -w-)P010 + jw' -u' )POl i + u'pm . 
Whenv>w>U ,P = - — -^^ . . . (lOe) 

As described above, according to this embodiment, 
since interpolation in the three-dimensional LUT can be 
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implemented by simple computations, high-speed 



processing can be easily realized. 



In the above description, a linear output is 



obtained. However, since the processing of this 
5 embodiment independently processes the respective 
output dimensions independently of the number of 
dimensions, this embodiment can be applied to 
conversion of an arbitrary output dimension. 



10 plurality of unit rectangular hexahedra, as shown in 
Fig. 6, and the shape of each unit rectangular 
hexahedron is determined by the grid position with 
respect to the respective dimensions. This embodiment 
can be directly applied to interpolation using an LUT 

15 in which the grid positions of a three-dimensional 
input are independently set, as shown in Fig. 6. 
However, if the same grid positions are set for all the 
dimensions, since the X-u' table can be used as the 
Y-v' and Z-w' tables, the processes for preparing these 

20 tables and a memory for storing them can be omitted. 

When this embodiment is applied to color conversion, if 
the input is one of RGB, CMY, and CIE XYZ of color 
spaces, the diagonal line of the LUT can be expressed 
by a gray line by setting grid positions of all the 

25 dimensions at identical positions, the precision of 

gray line interpolation can be improved. As described 



The actual three-dimensional LUT is composed of a 
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above, when the grid positions of all the dimensions 
are set at identical positions, since table preparation 
and a memory for storing the tables can be omitted, 
this embodiment can provide greater effects, 
5 [Arrangement of Conversion Apparatus] 

Fig. 8 is a block diagram showing the arrangement 
of the conversion apparatus. A CPU 102 controls the 
operation of the overall conversion apparatus via a bus 
107 in accordance with a program and data stored in a 

10 ROM 101, executes color conversion according to this 
embodiment for image data input via, e.g., an I/O 104 
or NIC (Network Interface Card) 105 using a RAM 103 as 
a work memory, and outputs the processing result via 
the I/O 104 or NIC 105 or stores it in a storage medium 

15 such as a hard disk (HD) 107. 

For example, to the I/O 104, an image I/O device 
such as a monitor (CRT, LCD, or the like) , printer, 
image reader, film reader, digital still caraera, 
digital video camera, or the like, or a storage device 

20 having a storage medium such as a magnetic disk, 
optical disk, or the like, is connected via a 
predetermined interface. Also, the apparatus can 
exchange image data with a computer to which such image 
I/O device or storage device is connected, via the NIC 

25 105. As such network, a network using Ethernet or FDDI 
(Fiber Distributed Data Interface) , and serial buses 
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specified by IEEE1394, USB (Universal Serial Bus), and 

the like are available. 

Fig. 9 is a flow chart showing conversion 

executed by the CPU 102. After a grid position is set 
5 in step SI, X-u', Y-v' , and Z-w' tables are prepared in 

steps S2 to S4. The prepared tables are stored in, 

e.g. , the RAM 103. 

Subsequently, image data is input in step S5, and 

u*, V*, and w' corresponding to the input image data 
10 are obtained using the prepared tables in step S6. The 

relationship among u', v', and w' is determined in step 

87, and converted image data is calculated using an 

equation corresponding to the determination result, 

i.e., one of equations (10a) to (lOf) in step S8. 
15 Steps S5 to S8 are repeated until all conversion 

processes of the image data are complete (not shown in 

Fig. 9) . , ■ ' 

Fig. 8 shows an example wherein the conversion 

apparatus of this embodiment is implemented by a 
20 computer such as a personal computer. However, this 

embodiment is not limited to such specific apparatus. 

For example, a DSP (Digital Signal Processor) may be 

provided in addition to the CPU 102, a program that 

implements the conversion shown in Fig, 9 may be 
25 supplied to the DSP, and the DSP may execute the 

conversion- Furthermore, the conversion apparatus of 
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this embodiment may be built in an image forming 
apparatus such as a printer, copying machine, or the 
like. In this case, X-u', Y-v', and Z-w' tables may be 
stored in an LUT RAM connected to an image data bus, 
5 and the conversion result may be computed by, e.g., a 
DSP on the basis of the outputs from the RAM, and may 
be output onto the image data bus. When the conversion 
process of this embodiment is applied to an image 
forming apparatus, that conversion process corresponds 
10 to various processes such as input masking, 

luminance-gray level conversion, UCR, output masking, 
gamma correction, gradation correction, color space 
conversion, color space compression/expansion, and the 
like - 

15 As described above, according to this embodiment, 

since a look-up table that obtains distance, which is 
normalized by a sufficiently large numerical value, 
from a grid point to an input value using the input 
value as an address is used, interpolation with 

20 practically high precision can be implemented by a 

computation volume equal to that required when all the 
grid spacings are equal to each other without limiting 
setups of grid positions. Hence, according to the 
conversion apparatus of this embodiment, flexible 

25 conversion characteristics can be implemented without 
increasing the computation volume and circuit scale. 
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This embodiment can especially reduce 
multiplications required for calculating weights when 
the input has two or more dimensions and the 
interpolation scheme is tetrahedral interpolation, thus 
5 effectively shortening the processing time and reducing 
the circuit scale. 

When a power of 2 is set as the sufficiently 
large numerical value for normalizing the distance 
between the input value and grid point, since divisions 
10 by all the weights done at the end of interpolation can 
be implemented by shift computations, thus further 
shortening the processing time and reducing the circuit 
scale . 



15 dimensions of the input are equal to each other, the 

processing steps and memory size required for 

implementing the present invention can be further 
reduced. 



20 effective to improve gray line interpolation precision, 
when the color space of the input is one of RGB, CMY, 
and CIE XYZ , 

To restate, according to the present invention, a 
data conversion method suitable for digital computation 
25 processing, a data conversion apparatus using the 

method, a recording medium used in the data conversion. 



Furthermore, when the grid positions for all the 



Note that this embodiment is particularly 
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and an image processing apparatus and method using the 
data conversion can be provided. 

As many apparently widely different embodiments 
of the present invention can be made without departing 
5 from the spirit and scope thereof, it is to be 

understood that the invention is not limited to the 
specific embodiments thereof except as defined in the 
appended claims. 
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